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Abstract 

The paper studies the excitability properties of a generalized FitzHugh-Nagumo model. 
The model differs from the purely competitive FitzHugh-Nagumo model in that it accounts for 
the effect of cooperative gating variables such as activation of calcium currents. Excitability 
is explored by unfolding a pitchfork bifurcation that is shown to organize five different types 
of excitability. In addition to the three classical types of neuronal excitability, two novel types 
are described and distinctly associated to the presence of cooperative variables. 

1 Introduction 

At the root of neuronal signaling, excitability is a dynamical property shared by all neurons, but 
its electrophysiological signature largely differs across neurons and experimental conditions. In the 
early days of experimental neurophysiology, Hodgkin [1] identified three distinct types of excitability 
(nowadays called Type I, II, and III) by stimulating crustacean nerves with constant current stimuli. 
The three types of excitability have long been associated to three distinct mathematical signatures 
in conductance-based models. They all can be described in planar models of the FitzHugh [2 type, 
which can be rigorously associated to the mathematical reduction of high-dimensional models (see 
for instance the planar reduction of Hodgkin-Huxley model in [3] and the excitability analysis in 0]). 
Reduced models have proven central to the understanding of excitability and closely related mecha- 
nisms such as bursting. However, understanding excitability in detailed conductance-based models 
remains a challenge, especially for neurons that exhibit transition between distinctively different 
firing types depending on environmental conditions. The present paper proposes a generalization 
of FitzHugh-Nagumo model that provides novel insights in the simple classification of excitability 
types. 

The proposed model is a "mirrored" version of FitzHugh-Nagumo model, motivated by the 
mathematical reduction of conductance-based models including calcium channels. A central obser- 
vation in [5] is that the cooperative nature of calcium channel activation makes their contribution to 
excitability fundamentally different from competitive gating variables, such as sodium inactivation 
or potassium activation. This chief difference is responsible for an alteration of the phase portrait 
that cannot be reproduced in FitzHugh-Nagumo model and that calls for a generalized model that 
motivates the present paper. 

We construct a highly degenerate pitchfork bifurcation (co-dimension 3) that is shown to orga- 
nize excitability in five different types. The three first types correspond to the types of excitability 
extensively studied in the literature. They are all competitive, in the sense that they only involve a 
region of the phase plane where our model is purely competitive. In addition, the model reveals two 
new types of excitability (Type IV and V) that match the distinct electrophysiological signatures of 
conductance-based models of high density calcium channels. We prove that these two new types of 
excitability cannot be observed in a purely competitive model, such as FitzHugh-Nagumo model. 

Global phase portraits of the proposed model are studied using singular perturbations theory, 
exploiting the timescale separation observed in all physiological recordings. An important result 
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of the analysis is that Type IV and V excitable models exhibit a bistable range that persists in 
the singular limit of the model. This, in sharp contrast to the Type I, II, and III. The result 
suggests the potential importance of Type IV and V excitability in bursting mechanisms associated 
to cooperative ion channels. 

The paper is organized as follows. The studied planar model, its geometrical properties, and 
the underlying physiological concepts are presented in Section [2j In Section [3j we construct the 
pitchfork bifurcation organizing the model and unfold it into the five types of excitability. Types 
I, II, and III are briefly reviewed in Section |4| Types IV and V are defined and characterized via 
phase plane and numerical bifurcation analysis in Section [5| By relying on geometrical singular 
perturbations, we provide in Section [6] a qualitative description of the phase plane of Type IV and 
V excitable models. Their peculiarities with respect to Type I, II, and III are particularly stressed. 
A summary and a discussion of our analysis are provided in Section [7J 



2 A mirrored FitzHugh-Nagumo model and its physiological 
interpretation 

The paper studies the excitability properties of the planar model 
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The model is reminiscent of the popular FitzHugh-Nagumo model of neuronal excitability: equation 
([iji) describes the fast dynamics of the membrane potential V, whereas ^p) describes the slow 
dynamics of the "recovery variable" n that aggregates the gating of various ionic channels. 

The voltage dynamics |lj,) are identical to the FitzHugh-Nagumo model, except that the 
quadratic term n 2 replaces the linear term n. This sole modification is central to the result of 
the present paper. The resulting nullcline V = "mirrors" along the V-axis the classical inverse 
V-shaped nullcline of FitzHugh-Nagumo model, as illustrated in Figure [TJ The figure illustrates 
the phase portrait of the model in the singular limit e = 0. The left and right phase-portraits are 
the unfolding of the transcritical bifurcation organizing center (see e.g. [51 Pages 104-105]) obtained 
for I app = I* := |. This particular value will help understanding excitability mechanisms at work 
in the situation I app > I* , illustrated in the right figure. 

Transcritical organizing center 
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Figure 1: V^-nullclines of |T]) for different values of I apP and e = 0. Stable fixed points are depicted as filled 
circles, unstable as circles, and bifurcations as half filled circles 

The recovery dynamics (JIJd) exhibit the familiar firts-order relaxation of ionic current to the 
static sigmoid curve illustrated in Figure [5J For quantitative purposes, the numerator in the right 
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hand side of ([2| can be picked larger or equal than 2 without changing the underlying qualitative 
analysis. The parameters (Vq, no) locate the relative position of the nullclines in the phase portrait. 
In particular, the region 

S„ := {(V,n) € M 2 : «£ (no, n + 2)} (3) 

is attractive and invariant for the dynamics 0. The parameter no slides up and down the "phys- 
iological window" of the recovery variable, whereas, Vq is the half-activation voltage. The voltage 
V* < V is defined as the voltage at which n 00 (F — Vq) has unitary slope. We adopt the conven- 
tional notation n for the recovery variable but will allow no < 0, which makes the range of n include 
negative values. This is purely for mathematical convenience and should not confuse the reader 
used to the physiological interpretation of n as a gating variable with range [0, 1]. 
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Figure 2: Dependence of the n- nullclines of Q on the parameters no and Vq. The location of the nullcline 
in the phase plane determines an attractive invariant region S no - At the voltage V*, the nullcline has 
unitary slope. 



Competitive and cooperative excitability 

In the terminology of [TJ, the variables V and n of model are competitive in the half plane n > 0, 
that is the model satisfies 

dV dn < 
dn dV 

whereas they are cooperative in the half plane n < 0, that is 

dV dn ^ 
dn dV 

A consequence of that observation is that model ([T]) is purely competitive when no > 0. The 
corresponding phase portrait is illustrated in Figure [3|\. It is reminiscent of FitzHugh-Nagumo 
model and will recover known types of excitability. In contrast, the model is neither competitive 
nor cooperative when no < 0. The corresponding phase portrait (Figure |3j3) is distinctly different 
from FitzHugh-Nagumo model and it will lead to the novel types of excitability studied in this 
paper. The role of the proposed mirrored FitzHugh-Nagumo model is to study the transition from 
a purely competitive model to a model that is neither cooperative nor competitive through a single 
parameter no. 

The physiology behind competitive and cooperative behaviors 

Competitive and cooperative behaviors model different types of (in)activation gating variables in 
conductance-based models. When the recovery variable is purely competitive (Figure[3]A_), it models 
the activation (resp. inactivation) of an outward (resp. inward) ionic current: a positive variation 
of V induces a positive variation of n and thus an increase of the total outward ionic current, 
i.e. n 2 , which is monotone increasing in S no . In contrast, when the recovery variable becomes 
cooperative (Figure [3^3), its role is reversed: it models the activation (resp. inactivation) of an 
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A. FHN type of excitability B. Novel type of excitability 




Cooperative Competitive Cooperative Competitive 




Figure 3: Phase portrait and outward ionic current n 2 for different positions of the n-nullcline and the 
associated invariant region S no . A) When S„ is fully contained in the half-plane n > the phase portrait 
is reminiscent of original FitzHugh-Nagumo model and the outward ionic current is monotone increasing. 
The recovery variable has a purely competitive role. B) When S no extends to the half-plane n < 0, the 
phase portrait exhibits new characteristics. The outward ionic current is not monotone, corresponding to 
two antisynergistic (competitive and cooperative) roles of the recovery variable. 

inward (resp. outward) ionic current: a positive variation of V induces a positive variation of n and 
thus a decrease of the total outward ionic current, which is now monotone decreasing. 

The seminal model of Hodgkin-Huxley only includes competitive slow gating variables: inacti- 
vation of sodium current and activation of potassium current. That is why the classical reduction 
of the Hodgkin-Huxley model leads to a FitzHugh-Nagumo type of phase portrait (Figure [3j\) . 
However, conductance-based models often include cooperative gating variables. An example of the 
latter is the activation of calcium currents included in many bursting conductance-based models 
(e.g. R15 neuron of Aplysia's abdominal ganglion 8 a , thalamo-cortical relay and reticular neurons 
[9lll0j. CA3 hippocampal pyramidal neuron [11)). Adding a calcium current in the Hodgkin-Huxley 
model is one natural way to obtain a reduced phase portrait as in Figure[3}3. This observation, firstly 
presented in [5] and reproduced in Figure |4j motivated the present study. 

3 A pitchfork bifurcation organizes different excitability types 

The model ([I]) has three free geometrical parameters (I a pp, no, Vo). The parameter is an ad- 
ditional parameter with respect to FitzHugh-Nagumo dynamics. The three parameters can be 
adjusted to create a codimension three bifurcation that will provide an organizing center for ex- 
citability. 

The degenerate bifurcation is illustrated in Figure [5] and is constructed as follows: 

1. The applied current is fixed at I app — I*, imposing the existence of the singularly perturbed 
transcritical bifurcation at the y-nullcline self- intersection (cf. Figure [l] center) . 

2. We fix n = «o(Vb) := — n-oo(— 1 — Vo) to force a nullcline intersection at the transcritical 
singularity at (—1,0). 

3. We fix V* = —1, or, equivalently, Vo = V * := — 1 + | log (—6 + \/35), so that the n-nullcline 
is tangent to the V-nullcline at the intersection. 

A local normal form of the model |T]) at this degenerate bifurcation is provided by the following 
lemma. 
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Figure 4: In original reduced Hodgkin-Huxley model a transcritical singularity exists in the non- 
physiological region of the phase plane. The addition of a calcium current makes this bifurcation phys- 
iological. 
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Figure 5: Phase portrait of Q for I app = I*, no = no(Vb), and V* = — 1. 



Lemma 1 There exists an affine change of coordinates that transforms |I]) into 



- (k(V )(v - U) + J ) 2 + Iapp - I* 



u = —eu + 0(v 2 ,vu, u )), 

where fc(Vb) := ^\p( — 1 — Vo) * s the slope of the activation function n^V — Vq) at V 
$o = n + n OQ (-l - V ). 

Proof The affine change of coordinates 

w := n — (5o, 
v = V + 1, 



(4a) 
(4b) 
-1, and 



transforms ([!]) into 



(W + So) 2 + Iapp ~ I* 



w = e(noo(v — 1 — Vb) — «oo(-l - Vo) - w), 



(5) 
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where So is defined as in the statement of the lemma. Because w\ v =w=o = 0, we extract the linear 
term in ^ and write 



v = — ^ - (id + S ) 2 + I app - I* 



a pi) 

2 „,„7, „r,2\ 



w = e(k(Vo)v — w + 0(v , vw, w )), 
Finally, the linear transformation u = v — k ^y o ^ transforms \Q\ into 



(6b) 



(k(V )(v - u) + So) 2 + I apP - I* 



u = —eu + 0(v 2 , vu, u 2 )), 

which proves the lemma. □ 
The dynamics Q have an exponentially attractive center manifold that is tangent at v = u = 
to center space {(y, u) : u — 0}. Ignoring higher order terms, the dynamics on the center manifold 
is given by 



v = (1 - k(V Q ) 2 )v 2 - — - 2k(V Q )vS - 5 2 + I app - I*. 



(7) 



Codimension 2 transcritical bifurcation for I app = I*, n = tiq(Vo), Vq 7^ V * 

Fixing I app — I* (Item 1.) and So — (Item 2.), but Vq 7^ V * , we obtain from ([7| the following 
relationships: 
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(8b) 



(8c) 



where ([8p) comes from the fact that, since Vq 7^ V *, k(Vo) 7^ 1, whereas (|8p) comes from the fact 
that k(Vo) 7^ for all Vq G K. Relation ^ are the defining conditions of a transcritical singularity 
(see e.g. Page 367]). For any value of the perturbation parameter So 7^ 0, there are two fixed 
points that exchange their stability at the bifurcation for So = 0. Figure [6] A,C illustrates this 
result. Note that we have to fix exactly two parameters (i.e. I app and 5q) to obtain ([8]), which 
confirms that bifurcation described by (l8|) has codimension 2. 



Codimension 3 pitchfork bifurcation for I app = I*, no = 71q(Vo), Vo = V * 
Adding the condition Vq = V * (Item 3.), we obtain the extra degeneracy condition 



d 2 i> 
dv 2 



„= 5 =o =2(l-fc(V ) 2 ) =0, 



and the associated bifurcation is in this case a codimension 3 pitchfork bifurcation (see e.g. O Page 
367]). In this case, positive perturbation of So leads to a unique stable fixed point, which split in 
three fixed points, the outer stable and the central unstable, at the bifurcation (see Figure [6] B). 



6 



A. lapp=l*,V >Vo* 




Figure 6: Phase portrait of |T]) for I apP = I*, different values of Vo, and with no as the bifurcation 
parameter. A,C: When Vo < Vq or Vo > V * as no decreases below no(Vb) a stable and an unstable 
fixed points exchange their stability in a codimension 2 transcritical singularity. B: When Vo = Vo* the 
bifurcation degenerates in codimension 3 pitchfork at which a stable fixed point splits in two stable fixed 
points (outer) and an unstable fixed point (inner). 



Unfolding in the plane I app = I* 

The parameter chart in Figure [7] unfolds the pitchfork bifurcation in the plane (Vo, no) for the fixed 
critical value I*. This bifurcation analysis reveals five qualitatively distinct regions denoted by I, 
II, III, IV, and V. The transition from Region I to Region II is through a saddle-node bifurcation 
at which the n-nullcline is tangent to the V-nullcline. See how the top right phase portrait is 
continuously deformed to the top left phase portrait in Figure [7J The transition from Region I 
to Region IV is through a transcritical bifurcation at which a saddle and a node exchange their 
stability. See how the top right phase portrait is continuously deformed to the bottom right phase 
portrait in Figure [7j A similar transition occur from Region IV to Region V. See how the bottom 
right phase portrait is continuously deformed to the bottom left phase portrait in Figure [7j Finally, 
the transition from Region V to Region II is through a saddle-node bifurcation. See how the bottom 
left phase portrait is continuously deformed to the top left phase portrait in Figure [7] 

Region III in Figure [7] is illustrated for future reference in the next section, but it should not 
be differentiated from Region II in the plane 1 = 1*, i.e. there is no bifurcation associated to the 
transition from Region II to Region III. 

The pitchfork bifurcation organizes excitability 

The relevance of the unfolding in Figure [7] for excitability is that the different regions correspond 
to different types of excitability for I app > I*. Regions I, II, and III correspond to Types I, II, and 
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Figure 7: Unfolding of the degenerate pitchfork bifurcation in the plane I ap p = !*■ Stable fixed points are 
depicted as filled circles, whereas unstable as circles. Saddle points are depicted as crosses. ■: pitchfork 
bifurcation. TC: transcritical bifurcation. SN: saddle-node bifurcation. 

Ill excitability identified in the early work of Hodgkin pQ and extensively studied in the literature 
since then. Those are the only types of excitability that can be associated to purely competitive 
models (i.e. no > 0) and they all have been studied in FitzHugh-Nagumo type phase portraits. 
They are briefly reviewed in Section [4] In contrast, Region IV and V correspond to new types of 
excitability that require the co-existence of competitive and cooperative ionic currents. They are 
studied in Sections [5] [6] 

Our analysis assumes a timescale separation e <C 1, reflecting the accepted strong separation 
between the fast voltage dynamics and sodium activation kinetics and the remaining slow gating 
kinetics. We focus on those bifurcations that persist in the singular limit e — > 0. 

Figure [§] summarizes the different types of excitability studied in the next section and their main 
electrophysiological signatures. 

4 Three types of excitability in competitive models 
Type I (SNIC) 

Fixing the parameters (Vo,no) in Region I of the parameter chart in Figure [7j the node and the 
saddle approach each other as the applied current I app > I* increases and eventually collide in a 
saddle-node bifurcation at the critical value I app = Isnic, as depicted in Figure [8] (top right). Un- 
der the timescale separation assumption (i.e. e <C 1), the center manifold of the bifurcation forms 
a homoclinic loop, as sketched in the figure. This bifurcation is commonly identified as saddle-node 
on invariant circle (SNIC). As I app is further increased, the fixed point disappears and the system 
generates a periodic train of action potentials. 

The excitability properties of model ([T]) near a SNIC bifurcation are commonly referred to as Type I 
excitability (see e.g. [TH Section 3.4.4], 0], and [TU Section 7.1.3] and references therein). The 
main properties associated to Type I excitability are as follows: 

All-or-none spike. The model has a well defined threshold (i.e. I app = Isnic) to generate 
action potentials. Above the threshold, the amplitude of the action potentials does not depend 
on the stimulus intensity. 

Low frequency spiking. The frequency of the limit cycle decreases to zero as I app \ Ig nic ■ 
From a computational point of view, this property permits to encode the stimulus intensity 
in the oscillation frequency. 
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Figure 8: Sketch of different types of excitable behaviors in the unfolding of the degenerate pitchfork 
bifurcation (■ in the central parameter chart) for I apP > I* and e< 1. For each type, the typical voltage 
time course and the phase portrait are sketched. The abbreviation SNIC denotes the saddle-node on 
invariant circle bifurcation. Stable fixed points are depicted as filled circles, unstable as circles, saddles 
as cross, and bifurcations as half-filled circles. The stable manifold of saddle points is depicted in green. 
The center manifold of the SNIC bifurcation is depicted in blue. The saddle homoclinic loop in Type IV is 
depicted in orange. 



9 



Examples of neurons exhibiting Type I excitability include: thalamocortical neurons with inacti- 
vated T-type calcium current (depolarized steady-state) [H] Figure 3], isolated axons from Carcinus 
maenas [TJ Class 1], regular spiking neurons in somatosensory cortex [T5], molluscan neurons [T^j . 

Type II (Hopf ) 

In Region II of the parameter chart in Figure [7j the nullclines intersect only once at a stable fixed 
point. As I app is increased, this fixed point looses stability in a Hopf bifurcation at I app = I Hopf, 
as illustrated in Figure [8] (top left). Above this critical input current, the system possesses a stable 
limit cycle surrounding the unstable fixed point. 

Excitability properties associated to a Hopf bifurcation are well known and define Type II excitabil- 
ity (see e.g. [TH Section 3.4.4], [3], and [T31 Section 7.1.3] and references therein). Fundamental 
properties of Type II excitable systems include: 

No threshold. When the bifurcation is supercritical, the amplitude of the limit cycle de- 
creases to zero as I app \ I Hopf-, which makes it hard to define a threshold for the generation 
of an action potential. Canard trajectories [17] are sometime considered as "soft" threshold 
manifold between small and large amplitude action potential. 

Subthreshold oscillations. When the applied current is slightly below the bifurcation 
values (i.e. I app < Inopf) the system trajectory relaxes to the fixed point with damped 
oscillations at the natural frequency of the Hopf bifurcation (i.e. the imaginary part of the 
eigenvalues at the bifurcation). 

No low frequency firing. The oscillation frequency is (almost) independent of the injected 
current and is equal to or larger than the natural frequency of the Hopf bifurcation. 

Frequency preference. Trains of small (< Inopf) amplitude inputs can induce spike if 
the intra stimulus frequency is resonant with the natural frequency of the Hopf bifurcation. 
This phenomenon is tightly linked to the presence of subthreshold oscillations and permits to 
detect the presence of resonant harmonics in the stimulus. 

Post-inhibitory spike. Transient negative current can induce an action potential in Type 
II excitable systems. 

Examples of neurons exhibiting Type II excitability include: isolated axons from Carcinus mae- 
nas [U Class 2] , fast spiking neurons in somatosensory cortex [15] , alpha moto- neurons [18] . 

Type III 

Type III excitability was only recently studied [19]. It can be thought as a less excitable variant 
of Type II excitability. In Region III of Figure [7| the half-activation voltage Vq is so negative that 
the stable focus never looses its stability as the applied current is increased. Nevertheless, the 
model is still excitable. For instance, as depicted in Figure [8] (bottom left), a positive current step 
instantaneously shifts the stable fixed point upright and an originally resting trajectory is attracted 
toward the right branch of the y-nullcline before relaxing back to rest. On the contrary, if the 
applied current varies slowly no action potential is generated. 

Particular neurocomputational properties of Type III excitability have recently been highlighted in 

Slope detection. Because a brutal variation of the applied current is necessary to excite the 
model, Type III excitable neurons acts as slope detectors with a high temporal precision 

Slope based stochastic resonance. In the presence of noise, Type III excitable models 
are most sensitive to the stimulus slope and frequency, rather than to its amplitude. The 
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associated stochastic resonance phenomenon (slope based stochastic resonance) exhibits dis- 
tinctly different filtering properties with respect to the classical stochastic resonance in Type 
I/II excitable models. 

Examples of neurons exhibiting Type III excitability include: squid giant axons (revised model) 
[50] , auditory brain stem [TH] , isolated axons from Carcinus maenas jT[ Class 3] . 

5 Two novel types of excitability in non- competitive models 
Type IV (singularly perturbed saddle-homoclinic) 

As illustrated in the bottom right phase portrait of Figure [7j Type IV excitability is the first 
excitability type that involves the "mirrored" shape of the voltage nullcline in ([lj: the stable 
node and the saddle lie on the lower cooperative branch of the V-nullcline. In particular, the 
hyperpolarized stable steady state of Type IV excitable models lies in the cooperative region of the 
phase portrait, i.e. where 

dVdn % 
dn dV 

As a consequence, excitability properties of this phase portrait cannot be studied in FitzHugh- 
Nagumo like models. 

Fixing the pair (Vo ; n o) hi Region IV of Figure [7j we obtain the bifurcation diagram illustrated 
in Figure [9] together with the associated phase portraits. The stable node looses stability in a 
saddle-node bifurcation at the critical value I app = Isn > I* ■ For I app > Isn, the model possesses 
a stable limit cycle that attracts all solutions (but the unstable focus). The spiking limit cycle 
disappears in a (singularly perturbed) saddle-homoclinic bifurcation at I app — Ish- The stable 
node attracts all solutions (but those on the stable manifold of the saddle) for I app < Ish- Based 
on geometrical singular perturbations, we provide in Section [6] a global phase portrait analysis of 
Type IV excitable systems. 

The chief properties of Type IV excitable systems are: 

Bistability. Type IV excitable models are bistable in the parameter range Ish < Iapp < Isn- 
a limit cycle attractor coexists with a stable fixed point. 

Spike- latency. When the applied current is abruptly increased slightly above Isn, the 
trajectory necessarily travels the narrow region between the nullclines. Since the vector field 
is small in that region, the first action potential is fired with a large latency. 

After-depolarization potential. After the limit cycle has disappeared in the homoclinic 
bifurcation, the trajectory converges to rest by following the attractive branches of the voltage 
nullcline, thus generating robust ADPs (see [5]). 

We stress that bistability, spike-latency, and ADPs are all direct consequences of the presence of a 
saddle point on the cooperative branch of the voltage nullcline: the basin of attraction of the stable 
node and the limit cycle are separated by the stable manifold of the saddle; spike-latency reveals 
the "ghost" of the center manifold of the saddle-node bifurcation; ADPs are generated along the 
hyperbolic invariant structure provided by the saddle stable and unstable manifolds (see Section [6] 
below) . 

Examples of neurons exhibiting Type IV excitability include: subthalamic nucleus neurons |21j . 
thalamo-cortical reticular and relay neurons with deinactivated T-type calcium current (hyperpo- 
larized state) [55J [53J , dopaminergic neurons [M[ [55] , superficial pyramidal neurons [55] . 

Type V (saddle-saddle) 

Type V excitability relates to Type IV as Type III does to Type II: similarly to Type IV, the 
hyperpolarized stable steady state of Type V excitable models lies in cooperative region of the 
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Figure 9: Bifurcations in Types IV excitable systems. Top: bifurcation diagram. Stable fixed point are 
drawn as solid lines, unstable as thin lines, and stable limit cycles as think lines. SH denotes the saddle- 
homoclinic bifurcation, SN the saddle-node bifurcation. Bottom: phase portraits. The stable manifold of 
saddle points is depicted in green. Unstable fixed point are depicted as circles, saddles as crosses, and stable 
fixed points as filled circles. Limit cycle attractors are depicted in blue. Sample trajectories are depicted 
as black oriented lines. 



phase plane. The distinct feature of the bottom left phase portrait of Figure [7] is the co-existence 
of two stable fixed points, a "down-state" and an "up-state" . The saddle stable manifold separates 
the two attractors. As I app is decreased below /*, the up-state eventually looses its stability in a 
saddle-node bifurcation at IsN,up, leaving the down-state globally asymptotically stable. Similarly, 
as I app is increased above J*, the down-state eventually disappears in a saddle-node bifurcation at 
I = IsN,down- But the up-state itself eventually looses its stability in a Hopf bifurcation at I app — 
I Hopf- The Hopf bifurcation can either take place beyond the bistable range [Isn, up > IsN, down ] j & 
situation illustrated in Figure 10 or it can take place within the bistable range [IsN,up, Isn, down], 



in which case, depending on I app , the stable down-state coexists with either a limit cycle attractor 
or a stable fixed point, a situation illustrated in Figure [TT] 

The main properties of Type V excitable models are similar to those of Type IV and can be 
summarized as follows: 

Bistability Type V excitable models are bistable in the range Isn,up < Iapp < Is N, down'- a 
stable down-state coexists with an up-state attractor that can be cither a stable fixed point 
or a stable limit cycle. 

Spike latency. Similarly to Type IV, the down-state looses stability in a saddle-node bifur- 
cation on the lower (cooperative) branch of V-nullcline, leading to a long latency before the 
convergence to the up-state. 

Plateau potentials. The up-state having a higher voltage with respect to the down-state, 
the transition between the two gives rise to plateau potentials either with or without spikes. 

Examples of neurons exhibiting Type V excitability include: olfactory bulb mitral cells [27] and 
striatal medium spiny neurons [25] . More examples are listed at Scholarpedia journal article [29]. 
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Figure 11: Bifurcations in Types V excitable systems. Legend as in Figure [9] except HB denoting the 
Hopf bifurcation. The same bifurcation is depicted in the phase portrait a half-filled circle. 
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6 Singularly perturbed global phase portrait analysis of Types 
IV and V excitable models 



Under timescale separation assumption (i.e. < e <C 1), Types IV and V excitable models exhibit 
four distinct signatures with respect to Types I and II: (i) the existence of a saddle- homoclinic 
bifurcation and (ii) ADPs in Type IV; (in) bistability and (iv) spike latency in both Types IV 
and V. A global phase portrait analysis based on geometrical singular perturbations provides an 
analytical explanation for the occurrence of these signatures solely in Types IV and V. 



6.1 Singularly perturbed saddle-homoclinic bifurcation and robust ADP 
generation in Type IV excitable models 

In this section we prove the existence of the saddle-homoclinic loop in Type IV excitable systems. 
We then rely on this analysis to provide a qualitative picture of the ADP generation mechanism in 
these models. In the remainder of the section we assume that the pair (Vo, lies in Region IV of 
Figure [7} 

We start by briefly recalling some basic results of geometrical singular perturbation theory, using 
([!]) as an explicit example. The interested reader will find in |30j an excellent introduction to the 
topic, and in [3TJ[32l|33] some recent extensions on which we rely for the forthcoming analysis. The 
time rescaling r :— et transforms ([!]) into the equivalent system 

eV = V -^--n 2 + I app (9a) 
n = e(noo(V -V ) + n -n), (9b) 

which describes the dynamics ([lj in the slow timescale r. In the limit e = 0, commonly referred 
to as the singular limit, one obtains from ([I]) and ^ two new dynamical systems: the reduced 
dynamics 

= V-^-^ + Iapp (10a) 
h = e( nQO (V- Vo)+n -n), (10b) 
which evolve in the slow timescale r, and the layer dynamics 

V = V-^--n 2 +I app (11a) 
h = 0, (lib) 



which evolve in the fast timescale t. Figure |l2[a) depicts the fast-slow dynamics (10 1,(11 1. The 
main idea behind geometrical singular perturbation theory is to combine the analysis of the reduced 
and layer dynamics to derive conclusions about the behavior of the nominal system, i.e. with e > 0. 



The reduced dynamics ( 10 1 is a dynamical system on the set 



S := j(V» G R 2 : V - ^- - n 2 + I app = 

usu ally called the critical manifold. The points in Sq are indeed critical points of the layer dynamics 
( j 1 1 [ ) - More precisely, portions of So on which §^ is non- vanishing are normally hyperbolic invariant 
manifolds of equilibria of the layer dynamics, whose stability is determined by the sign of 

Conversely, points in So where ^ = constitute degenerate equilibria. In particular, the layer 
dynamics (11) exhibits, for I app — I* , two degenerate equilibri^ As depicted in Figure 12 'a) they 

1 A third degenerate equilibria is at the minimum of the V-nullcline that is specular to the maximum with respect 
to the line v = 0, but it plays no dynamical role and is not considered here. 
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(c) (d) 



Figure 12: Phase-portrait of ([T]) for (no,Vo) in Region IV of the parameter chart in Figure pf| 
and different values of e and /. (a): fast-slow dynamics of Q for e = and 1 = 1*. The attractive 
(resp. repelling) branches of the critical manifold So above/below the transcritical singularity are denoted 
by Sa / Sa (resp. S^ / S^). (b): Continuation of the slow attractive S a and repelling S* manifolds for e > 
and I apP < I* + I c (^/e), where 7 c (y / e) is defined as in Theorem[l] (c): Continuation of S% and 5y for e > 
and I = I* + /oCv^)- (d): Continuation of Sa and S* for e > and I > I* + I c {^/e). 
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are given by the self-intersection of the l^-nullcline, which is the transcritical organizing center in 
Figure [l] and by the fold singularity at the maximum of the upper branch of the y-nullcline. 

The basic result of geometrical singular perturbation theory, due to Fennichel [33], is that, for 
e sufficiently small, non-degenerate portions of So persist as nearby normally hyperbolic locally 
invariant manifolds S € of ([!]). More precisely, the slow manifold S e lies in a neighborhood of So of 



radius 0(e)- The dynamics on S e is a small perturbation of the reduced dynamics (10). We point 
out that S £ may not be unique, but is determined only up to 0(e~ c / e ), for some c > 0. That is, two 
different choices of S e are exponentially close (in e) one to the other. Since the presented results 
are independent of the particular S e considered, we let this choice be arbitrary. The trajectories of 
the layer dynamics perturb to a stable and an unstable invariant foliations with basis S £ . 

The analysis near degenerate points is more delicate. Only recently some works have treated 



this problem in its full generality for different types of singularities [3T1 [32J [33] . Figure 12 (b),(c),(d) 
sketch the extension of the attractive slow manifold 5* after the fold point, and the three possible 
ways in which and the repelling slow manifold S* can continue after the transcritical singularity, 
depending on the injected current. 



The result depicted in Figure 12 relies on the following analysis, adapted from 



Let A := {(V,n) el 2 : V- < V < V + , n = p}, be the section depicted in Figure 12 where 
p < and \p\ is sufficiently small, and V-,V+ are such that A n S~ ^ 0. For a given e > 0, 
let q a ^ t := AflS^ and q r ^ :— AnSJ be the intersections, whenever they exist, of respectively 
the attractive and repelling invariant submanifolds and S% with the section A. The following 
theorem reformulates in a compact way the discussion contained in Remark 2.2 and Section 3 of 
^lfor the system with inputs In). 



Theorem 1 (Adapted from |32j) Consider the system Q). Then there exists e > and a 
smooth function I c {^/e), defined on [0, eo] and satisfying I c (0) = 0, such that, for all e G (0, eo], the 
following assertions hold 

1- Qa.e = q r ,e if and only if I app = I* + I C (V*) 

2. there exists an open interval A 3 I* + / c (y / e), such that, for all I € A, it holds that AnS^ ^ 0, 
AnS e r ^ 0, and 

d 

(q a ,e ~ q r ,e) > 0. 



Figure fl2] illustrates this result. 



Remark 1 The function I c {s/e) is related to the function A c (y / e) defined in [32j Remark 2.2] by 
^c(v^) := eA c (v / e). Similarly, given e > 0, the parameter I app appearing in Theorem [l] is just the 
re-scaling I app = eX + I* of the parameter A appearing in [32, Remark 2.2 and Sections 3]. 

Theorem [T] implies the existence of the saddle- homoclinic bifurcation in the mirrored FitzHugh- 
Nagumo model with parameters (Vo,no) belonging to Region IV of Figure [j] In this case, as 
illustrated in Figure [l2|[b,c,d), the slow attractive manifold (resp. slow repelling manifold Sf.) 
coincides with the unstable manifold W u (resp. stable manifold W s ) of the saddle point, as it can 
be proved via simple qualitative arguments. Thus, for I app < I* + 7 c (i/e), the unstable manifold 
W u continues after the transcritical singularity on the left of W s , toward the stable node. See 
Figure [l2^b) . For I app = I* + 7 c (y / e), W u extends after the transcritical point to W s , forming the 
saddle-homoclinic trajectory, as sketched in Figure [8] and Figure [l2)[c) . For I app > I* + I c (^/e), the 
unstable manifold of the saddle W u continues after the transcritical singularity on the right of W s , 
and spirals toward an exponentially stable limit cycle, whose existence can be proved with similar 
geometrical singular perturbation theory arguments (see for instance |31j). This situation is the 



one depicted in Figure 12 [d). 

The existence of the saddle-homoclinic bifurcation near the singular limit e = is a direct 
consequence of the existence of a singular connection between the stable and unstable manifold of 



The first author is thankful to Prof. Szmolyan for his useful comments. 
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Figure 13: Singularly perturbed non-trivial saddle-homoclinic connection near the transcritical singularity. 
This loop persists in the non-singular limit. 



the saddle, as sketched in Figure 13 The existence of the saddle-homoclinic bifurcation can be seen 
as the persistence of this singular saddle-homoclinic loop in the non-singular dynamics. The same 
loop being absent in neuron models with an A-shaped nullcline, as discussed below, we identify 
with the singularly perturbed saddle-homoclinic bifurcation sketched in Figure [8] (bottom right) as 
the signature of Type IV excitable models. 

When I app < I* + I c (yfe) the above analysis ensures the robust generation of ADPs. The 
trajectory relaxation to the stable fixed point is indeed guided by the normally hyperbolic attractive 
manifold S^. As illustrated in Figure 12 'b), the voltage is not monotone as the trajectory slides 
along this manifold, corresponding to an ADP. The normally hyperbolicity of ensures that the 
ADP generation is robust to external perturbations. Some consequences for the modeling of neurons 
exhibiting ADPs can be found in [5]. 



Absence of singularly perturbed saddle-homoclinic bifurcation and ADPs in competi- 
tive models 

When the stable fixed point and the saddle belong to the upper A-shaped competitive branch of 
the voltage nullcline (i.e. Type I excitability), the model cannot exhibit saddle-homoclinic bifurca- 
tion. As sketched in Figure |14| in the fast-slow dynamics there are no nontrivial saddle-homoclinic 
connections. By persistence arguments 34. 31 , this implies that the nontrivial intersection of the 
stable and unstable manifolds of the saddle is empty also away from the singular limit. 

The prediction of singular perturbation theory of course does not contradict the existence of a 
saddle-homoclinic bifurcation in competitive models, it only precludes it for sufficiently small values 
of e. A well-known saddle-homoclinic bifurcation is described in Morris-Lecar model [35] for the 
barnacle giant muscle fiber. The model is purely competitive but a saddle-homoclinic bifurcation 
is possible around the N-shaped voltage nullcline (see for instance (T2], Section 3.4.3]). However, 
such a bifurcation cannot persist with a strong timescale separation, which suggests that it might 
be less relevant in the context of neuronal modeling. 




Figure 14: Absence of a non-trivial singularly perturbed saddle-homoclinic connection around the N- 
shaped V-nullcline branch. Such connection is absent also for sufficiently small e > 0. 
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Figure 15: Existence of a singularly perturbed saddle separatrix W a on the lower branch of the voltage 
nullcline in Type IV and V excitable systems for I > I* and different nullcline intersections. Such object 
persists in the non-singular limit as the saddle stable manifold Wl which separates the basin of attraction 
of the lower and upper attractors. 



The generation of ADPs must also be excluded when the stable fixed point belongs to the upper 
branch of the voltage nullcline. The relaxation to rest is indeed guided by the left attractive branch 
of the voltage nullcline, along which the voltage is monotone. ADPs can be generated within purely 
competitive models only by resorting to a non-physiological state reset mechanism JT3J Section 8.3]. 



6.2 Singularly perturbed bistability in Type IV and V excitable systems 

The geometrical singular perturbation machinery introduced above can be used to prove the per- 
sistence of bistability near the singular limit e = for both Type IV and V excitable systems. 
The analysis is sketched in Figure [15} The underlying key ingredient is the existence of singularly 
perturbed separatrix W% passing by the saddle point. This object persists in the non-singular limit 
as a normally hyperbolic saddle stable manifold that separates the stable down-state from the 
up-state attractor (see the bottom center plot in Figures [9 10 and 11). The phase portraits in the 
non-singular limit can all be deduced by the results in [32J [33] . 



Absence of singularly perturbed bistability in competitive models 

As opposed to Types IV and V, bistability is not observed in Types I, II, and III excitable models 
near the singular limit e = 0. For Type I this can be deduced by the absence of a singularly 
perturbed separatrix (cf. Figure 14 1 and standard persistence arguments. For Type II, one can 
invoke the fact that the bistable range of the subcritical Hopf bifurcation shrinks to zero in the 
singular limit (see for instance [33l Theorem 3.1]). Type III excitable models are, by definition, 
always globally asymptotically stable. 



6.3 Singularly perturbed spike latency in Type IV and V excitable sys- 
tems 

Spike latency appears when the trajectory travels the "ghost" of the center manifold of the saddle 
node bifurcation. It is prominent in Type IV and V excitable models, since in these models this 
center manifold is attractive. Indeed, the recovery variable nullcline being strictly monotone in- 
creasing, the saddle-node bifurcation lies on the lower attractive branch of the voltage nullcline (see 
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Figure 16: Mechanism of spike latency in Types IV and V excitable models. Left: slow-fast dynamics near 
the saddle-node bifurcation. Center: stable Wl and center Wl manifolds of the bifurcation for e > 0. Right: 
the ghost of Wl (depicted as the blue rectangle) is attractive for I > Isn corresponding to prominent and 
robust spike latency. 




s=0,I = I SN e>0,I = I SN 8=0,I>I SN 

Figure 17: Absence of spike latency in Type I excitable models. Left: the saddle-node bifurcation in Type 
I excitable models with its center Wl and stable Wl manifolds. Center: the stable and the saddle fixed 
points before the bifurcation. Right: the ghost of Wl (blue rectangle) after the biurcation and the behavior 
of an originally resting trajectory. 



Figure 16). By standard persistence arguments, this implies that its center manifold is strongly 
attractive near the singular limit. A consequence of this attractiveness is that, after the bifurcation 
at I app — Isn, an originally resting trajectory is attracted toward the ghost of Wl between the two 
nullclines, where the vector field magnitude is proportional to I app — Isn- The passage time thus 
diverges to infinity as I app \ Isn, corresponding to a prominent and robust spike-latency. 

Absence of singularly perturbed spike-latency in competitive models 

As opposed to Types IV and V, the saddle node bifurcation of Type I excitable models lies on 
the upper repulsive (right) branch of the voltage nullclinc (recall again that the recovery variable 
nullcline is strictly monotone increasing). Hence, attractiveness of the bifurcation center manifold, 
and thus of its ghost, is lost near the singular limit. Even though for finite e this center manifold can 



be attractive (see Figure 17 1, its ghost (between the two nullclines) does not attract an originally 
resting trajectory: the trajectory necessarily travels below the recovery variable nullcline where the 
vector field has finite magnitude independently of the distance from the bifurcation. These simple 
arguments show that competitive excitable models and, in particular, Type I excitable models can 
not exhibit spike latency. 



7 Discussion 

Mirroring the FitzHugh-Nagumo equation accounts for cooperative gating 
variables 

Mirroring the FitzHugh-Nagumo equation was motivated by the inclusion, in a simple model of 
neuronal excitability, of the transcritical bifurcation recently observed in a planar reduction of 
the Hodgkin-Huxley model augmented with an activating calcium current. This purely geometrical 
construction is actually tightly linked with the underlying electrophysiology: in the upper FitzHugh- 



19 



Nagumo like part of the phase portrait the model is competitive, as it is the original Hodgkin-Huxley 
model; in the mirrored part it is cooperative, accounting for cooperative gating variables, such as 
the activation of calcium current. To the best of our knowledge, this unified picture is not present 
in existing planar models of neuronal excitability, that are purely competitive. 



Cooperativity unravels a pitchfork bifurcation organizing old and new 
types of excitability 

The distinctive effects of cooperative gating variables on neuronal excitability are widely studied 
in high dimensional conductance-based models and in in vitro recordings. However, these same 
signatures can not be reproduced in competitive models and the underlying dynamical mechanisms 
have remained obscure to date. 

The unfolding of a pitchfork bifurcation organizing the proposed model reveals how the inclusion 
of cooperative variables changes neuronal excitability. The obtained parameter chart recovers the 
three known types of (competitive) excitability and unmasks two novel types that we naturally 
defined as Types IV and V. The defining condition of Types IV and V excitability is cooperativity, 
z c 

dVdn % 
dn dV 

at the hyperpolarized stable steady state. This sole condition ensures the presence of the electro- 
physiological signatures of cooperative gating variables, in particular, bistability. In a forthcoming 
publication, we will show that this defining condition also holds in higher-dimensional models and 
is not an artifact of planar reductions. 



Bistability and bursting in singularly perturbed excitable models 

A distinct feature of Types IV and V excitable models with respect to the three other types is the 
existence of a finite range of bistability. Bistability has been described in the context of Types I 



and II excitability as well (see e.g. [1]), but it was shown in Section 6.2 that, in all these situations, 
the bistability range shrinks to zero as the timescale separation is increased. In contrast, the stable 
manifold of the saddle point that separates the two basin of attraction in Types IV and V excitability 
is a robust (hyperbolic) geometric object that persists in the singular limit. Because neurons do 
exhibit a pronounced timescale separation, the robustness of the bistable range in Types IV and V 
is thought to be an important feature of excitable models that are not purely competitive. 

The relevance of bistability in excitable models lies in its relevance for model bursting. Bursting 
is typically the result of a slow adaptation variable that modulates the applied current across 
the bistability range, creating a hysteresis loop between the stable down-state and the up-state 
attractor. An important conjecture derived from our analysis is that bursting will persist near the 
singularly perturbed limit of model (JlJ only in Type IV and V excitability, that is only in the 
presence of cooperative ionic channels. In other words, Types IV and V excitability would be the 
essential sources of bursting in singularly perturbed models. 



Methods 

The parameter chart in Figures [7] and [8] have been numerically drawn using MATLAB^and modified 
with the Open Source vector graphics editor Inkscap^] 
Phase portrait were hand-drawn with Inkscape. 

The bifurcation diagrams in Figures [9j [ToJ and [IT] have been obtained with XPP environment^] 



"http : //www .mathworks . com 

4 http : //inkscape . org 

^http : //www . math . pitt ■ edu/~bard/xpp/xpp . html 
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